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NONLINEAR  TEMPERATURE  DEPENDENT  FAILURE  ANALYSIS 
OF  FINITE  WIDTH  COMPOSITE  LAMINATES 

ABSTRACT 

A  quasi-three  dimensional,  nonlinear  elastic  finite  element  stress 
analysis  of  finite  width  composite  lasiinates  including  curing  stresses 
is  presented. 

Cross-ply,  angle-ply,  and  two  quasi-isotropic  graphite/epoxy 
laminates  are  studied.  Curing  stresses  are  calculated  using  temperature 
dependent  elastic  properties  that  are  input  as  percent  retention  curves, 
and  stresses  due  to  mechanical  loading  in  the  form  of  an  axial  strain 
are  calculated  using  tangent  moduli!  obtained  by  Ramberg-Osgood 
parameters.  It  is  shown  that  oj(ring  stresses  and  stresses  due  to 
tensile  loading  are  significant  as  edge  effects  in  all  types  of 
laminate  studies. 

The  tensor  polynomial  failure  criterion  is  used  to  predict  the 
initiation  of  failure.  The  mode  of  failure  is  predicted  by  examining 
individual  stress  contributions  to  the  tensor  polynomial.  Failure 
is  predicted  to  always  initiate  at  the  free  edge,  but  not  always  at 
ply  interfaces.  The  location  and  mode  of  failure  is  shown  to  be 
laminate  dependent. 
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1 .  INTROOUCTION 


The  fabrication  of  laminated  fibrous  composites  involves  curing 
the  fiber-matrix  system  at  an  elevated  temperature.  Cure  temperatures 
for  resin  matrix  systems  vary  from  350°F  for  epoxies  to  650°F  for 
polyimides.  The  mismatch  of  thermal  expansion  coefficients  between 
fiber  and  matrix  (or,  alternatively,  the  orthotropic  properties  of  the 
lamina)  coupled  with  the  large  temperature  drop  from  the  maximum  cure 
temperature  can  result  in  relatively  high  residual  curing  stresses  in 
laminates  at  room  temperature.  These  thermal  stresses  are  often  large 
enough  to  cause  transverse  microcracking  or  ply  delamination  prior  to 
the  application  of  mechanical  load.  Residual  thermal  stresses  are  also 
present  in  metal  matrix  composites  such  as  boron-aluminum,  in  which 
their  effect  is  manifested  as  yielding  of  the  matrix  material. 

The  development  of  residual  stresses  in  composites  does  not  have 
a  direct  counterpart  in  homogeneous,  isotropic  media  and  there  are 
relatively  few  studies  of  the  subject  reported  in  the  literature.  All 
the  studies  are  based  upon  the  assumption  that  the  total  strain  is  the 
sum  of  two  distinct  parts:  the  mechanical  strain  which  is  related  to 
the  stresses  through  the  constitutive  equation,  and  the  "free"  thermal 
strain  which,  in  itself,  does  not  cause  stress  in  the  laminate. 

The  purpose  of  this  study  is  to  analyze  the  thermo-mechanical 
response  of  resin  matrix  composites  and  to  predict  the  occurrence  and 
mode  of  first  failure  in  finite  width  laminates.  Previous  researchers 
have  proposed  various  methods  for  calculating  residual  stresses,  and 
there  have  been  a  few  studies  of  stress-strain  response  to  mechanical 
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load  which  included  residual  stresses,.  Most  previous  studies  typically 
perform  the  residual  thermal  analysis  and  the  mechanical  load  analysis 
separately  assuming  linear  elastic  behavior.  The  principal  of  super¬ 
position  is  used  to  predict  the  combined  effect  of  mechanical  load  and 
curing  stresses.  The  present  study  treats  the  thermal  and  mechanical 
behavior  separately,  but  does  not  make  the  assumption  of  linear 
elasticity.  The  residual  stress  field' therefore  cannot  be  superposed 
on  the  mechanical  load,  but  is  used  as  an  initial  condition.  Special 
attention  is  given  to  the  influence  of  edge  effects  on  the  stress  field 
and  the  occurrence  of  first  failure. 

The  finite  element  program  N0NC0M2  [1,2,3]  was  modified  for  this 
analysis,  The  efficiency  and  capability  of  the  program  were  increased 
so  as  to  handle  a  larger  number  of  nodes  (with  a  choice  for  an  in-core 
or  an  out-of-core  equation  solver)  and  a  detailed  failure  analysis 
using  the  tensor  polynomial  failure  criterion  for  predicting  first 
failure.  The  modified  program  is  designated  N0NC0M3.  Results  were 
obtained  for  cross-ply,  angle-ply  and  quasi-isotropic  laminates  of  T300/ 
5208  graphite  epoxy. 


2. 


LITERATURE  REVIEW 


Most  previous  theoretical  studies  are  lamination  theory  solutions. 
They  are  based  upon  the  classical  plate  theory  assumptions  and,  there¬ 
fore,  valid  only  in  interior  regions  away  from  the  free  edges  of  finite 
width  laminates;  they  yield  only  laminate  stresses.  However,  failure 
in  laminates  is  often  observed  to  initiate  at  the  free  edges  [4,5]  and, 
therefore,  the  stress  distribution  there  is  of  paramount  interest. 

The  importance  of  edge  effect  during  thermal  loading  of  graphite- 
polyimide  was  clearly  demonstrated  by  Herakovich,  Mills  and  Davis  [5]. 

Tsai  [6]  presented  a  thermoelastic  formulation  for  calculating  lam¬ 
inate  thermal  stresses  in  1955.  This  study  presents  the  basic  lamination 
theory  development  for  calculating  residual  stresses-.  A  micromechanical 
procedure  for  calculating  residual  thermal  stresses  was  outlined  by 
Hashin  [7].  One  of  the  earlier  reported  analytical  predictions  of 
residual  stresses  using  lamination  theory  is  a  study  by  Chamis  [8], 
in  which  he  analyzed  laminates  of  different  material  systems,  stacking 
sequences  and  fiber  volume  fractions.  Extensive  experimental  studies 
were  conducted  at  the  IIT  Research  Center  by  Daniel  and  Liber  [9]. 

They  reported  thermal  stresses  based  upon  measured  strains  and  tempera¬ 
ture  dependent  constitutive  relations. 

Herakovich  [10]  was  apparently  the  first  to  consider  thermal 
edge  effects  in  composites.  He  analyzed  laminates  of  boron-epoxy  and 
aluminum  using  finite  elements,  and  considered  stress  distributions 
due  to  thermal  and  mechanical  loads.  The  analysis  included  interlaminar 
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stresses  but  was  linear  elastic  with  constant  room  temperature  properties. 
A  nonlinear  elastic  finite  element  analysis,  which  included  thermal 
effects  and  temperature  dependent  properties  was  conducted  by  Renieri 
and  Herakovich  [1],  but  residual  stress  predictions  formed  only  a 
limited  part  of  the  study  and  the  finite  element  mesh  used  was  quite 
coarse.  Their  basic  formulation  will  be  used  in  the  present  analysis 
with  a  much  finer  mesh,  an  improved  equation  solver  and  a  failure 
criterion  added  to  the  analysis  capability. 

Hahn  and  Pagano  [11]  pointed  out  the  necessity  for  the  inclusion  of 
terms  corresponding  to  the  stress  and  temperature  dependence  of  proper¬ 
ties.  They  developed  a  'total  strain'  theory,  in  which  the  strains  and 
stresses  are  calculated  using  temperature  dependent  elastic  properties 
at  the  temperature  of  interest. 

Daniel,  Liber  and  Chamis  [12]  developed  a  technique  to  measure 
residual  strain  by  embedding  strain  gages  between  plies  in  laminates. 

They  used  this  technique  for  measuring  curing  strains  in  boron/epoxy 
and  S  glass/epoxy,  and  calculated  stresses  using  temperature  dependent 
constitutive  relations.  Thermal  cycling  suggested  that  residual 
stresses  during  the  curing  process  were  primarily  due  to  thermal 
mismatch  between  adjacent  plies. 

Chamis  and  Sullivan  [13]  outlined  a  procedure  for  nonlinear  analysis 
of  laminates  with  residual  thermal  stresses.  The  laminate  was  loaded  in 
increments,  using  stresses  calculated  in  the  most  recent  load  step  to 
calculate  elastic  moduli  for  the  next  load  step.  Micromechanics  was 
used  to  predict  lamina  properties,  which  were  used  in  the  lamination 
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theory  analysis. 

Hahn  [14]  concluded  that  the  stress  free  temperature  in  laminates 
is  lower  than  the  cure  temperature.  The  method  outlined  in  [11]  was 
used  to  calculate  residual  strains  which  were  compared  to  experimentally 
determined  strains.  Daniel  and  Liber  [15]  investigated  the  effect  of 
stacking  sequence  on  residual  stresses  in  graphite/polyimide  laminates. 
The  strains  were  determined  experimentally,  and  the  stresses,  calculated 
using  constitutive  relations,  were  found  to  be  close  to  the  transverse 
strength. 

Wang  and  Crossman  [16]  studied  edge  effects  due  to  thermal  loading 
on  some  specific  laminates.  They  predict  a  peculiar  behavior  for  a 
[±45]^  laminate,  with  the  existence  of  'stiff  tensile  and  'soft' 
compressive  zones  in  the  laminate. 

A  report  by  Chamis  [17]  summarized  work  done  at  the  NASA  Lewis 
Research  Center  on  angle  ply  laminates  ove»“  a  period  of  eight  years. 

The  effect  of  curing  stresses  on  laminate. warpage  and  fracture  was 
studied  experimentally  and  analytically  using  lamination  theory. 

Pagano  and  Hahn  [18]  used  the  procedure  described  in  [11]  to  calcu¬ 
late  residual  thermal  stresses,  and  studied  their  effect  on  failure  in 
laminates.  The  curing  stresses  were  found  to  influence  first  failure  in 
laminates  greatly,  often  reducing  the  applied  load  to  failure  by  about 
half.  They  note  that  the  interlaminar  normal  stress  is  significant 
in  some  stacking  sequences,  especially  at  free  edges,  and  that  this 
would  result  in  failure  initiating  at  loads  much  less  than  their  calcu¬ 
lated  values.  Their  analysis  is  based  on  lamination  theory  and  thus 
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does  not  treat  the  free  edge  problem  in  any  detail. 

Farley  and  Herakovich  [19],  using  a  finite  element  analysis,  compared 

boundary  layer  stress  distributions  due  to  mechanical,  thermal,  and 
moisture  loads  in  finite  width  laminates.  Each  type  of  load  was  analy¬ 
zed  separately;  the  study  concentrated  bn  the  response  of  laminates  to 

different  moisture  gradients  in  the  boundary  laye  . 

Kim  and  Hahn  [20]  published  results  of  acoustic  emissions  of 
laminates  subject  to  -.ochanical  loads.  Curing  stresses  were  included 
in  the  lamination  theory  development  for  predicting  stress  at  which 
first  failure  occurred. 


3.  THEORETICAL  BACKGROUND 


The  problem  under  consideration  is  the  stress  analysis  of  synmetric 
laminates,  including  thermal  and  free  edge  effects.  In  this  study  the 
nonlinear  analysis  for  both  mechanical  and  thermal  loading  is  performed 
using  an  incremental  procedure.  The  loading  is  approximated  by  a  finite 
number  of  load  steps  and  each  step  is  treated  as  a  linear  problem.  The 
applied  load,  whether  mechanical  or  thermal,  is  assumed  to  be  steady 
and  uniform  across  the  laminate. 

3.1  General  Formulation 

A  typical  balanced,  symmetric  laminate  is  shown  in  Fig.  1.  The 
behavior  of  the  laminate  can  be  assumed  to  be  independent  of  the  x 
coordinate  if  b  and  H  are  small  compared  to  L.  As  shown  by  Pipes  and 
Pagano  [21]  the  linear  strain  displacement  relations  can  be  integrated 
and  manipulated  to  yield  the  following  displacement  field  over  the 
cross-section  of  the  laminate. 

u  -  — (CiZ+C2)y  ( G^y"*"CgZ+Cg )x  +  U(y,z) 

V  =  (C^z+C2)x  -  C4  V(y,z)  (3.1) 

x2 

w  =  -C^xy  +  CyX  -  Cg  ^  +  Cg  +  W(y,z) 

The  displacement  field  has  the  following  symmetry:  with  respect  to 
the  x-y  plane, 
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u(x,y,z)  =  u(x,y,-z) 

v(x,y,z)  =  v(x,y,-z)  (3.2a) 

w(x,y,z)  =  -w(x,y,-z) 
with  respect  to  the  x-z  plane, 

v(x,y,z)  =  -v(x,-y,z) 


(3.2b) 


w(x,y,z)  =  w(x,-y,z) 

It  has  been  experimentally  observed  [22]  that  at  z=±H 

u(x,y,±H)  =  -u(x,-y,±H)  (3.3) 

As  the  thickness  of  the  laminate  is  small,  it  can  be  assumed  that 

u(x,y,z)  =  -u(x,-y,z)  ,  (3.4) 

These  symmetries  simplify  the  displacement  field  to 

u  =  CgX  +  U(y,z) 


V  =  V(y,z)  (3.5) 

w  =  W(y,z) 

The  analysis  can  now  be  restricted  to  one  quarter  of  the  cross- 
section  (Fig.  2)  with  the  following  boundary  displacement  constraints: 


V(0,z)  =  0 
W(y,0)  =  0 


(3.6) 


The  following  stress  boundary  conditions  complete  the  boundary  value 
problem. 


FIGURE  2.  BOUNDARY  CONDITIONS  ON  THE  QUARTER  SEC  1  ION  OF 
THE  LAMINATE 


n 


=  T2^(x,y,±H)  =  a^Cx.y.+H)  =  0 

(3.7) 

Tj^y{x,±b,z)  =  T^y(x,±b,2)  =  Oy.'x.ib.z)  =  0 

The  individual  laminae  are  orthotropic,  having  a  stress  strain 
relation  with  9  independent  constants.  When  referred  to  the  laminate 
axis,  the  stress  strain  relation  transforms  to  (Appendix  A) 


(3.8) 


3.1.1  Finite  Element  Formulation 

This  boundary  value  problem  is  cast  in  the  finite  element  frame¬ 
work.  The  cross  section  is  subdivided  into  triangular  elements,  and  the 
displacement  field  is  assumed  to  vary  linearly  within  each  element.  The 
Elemental  displacement  field  is  represented  in  terms  of  the  nodal  coor¬ 
dinates  and  the  nodal  displacements.  The  total  potential  energy,  con¬ 
sisting  of  the  strain  energy  and  the  potential  of  external  forces,  is 
written  for  each  element  in  terms  of  the  nodal  displacements  and  forces. 
The  potential  energy  is  then  minimized  with  respect  to  the  nodal 
displacements  to  obtain  a  linear  set  of  equations  relating  nodal 
displacements  to  nodal  forces  through  the  element  "stiffness  matrix". 
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These  elemental  stiffness  matrices  are  assembled  to  form  a  system  of 
equations  in  the  unknown  nodal  displacements.  The  system  of  equations 
is  solved  after  imposition  of  boundary  conditions.  The  strains  and 
stresses  in  each  element  are  calculated  from  the  displacements  of  the 
element  nodes,  the  strain  displacement  relations  and  the  constitutive 
equations. 

3.2  Mechanical  Loading 

Let  the  laminate  in  Fig.  1  be  loaded  with  a  uniform  strain  5  in 
the  X  direction.  The  displacement  field  over  an  element  at  a  cross 
section  x=Xi  becomes 

u  =  a^  +  a2y  +  a3Z  +  cx^ 

V  =  a^  +  a^y  +  a^z  (3.9) 

w  =  a?  +  agy  +  agZ 

When  the  parameters  a-j-ag  are  functions  of  the  nodal  coordinates  and 
displacements.  As  the  laminate  behavior  is  independent  of  the  x 
coordinate,  x^  is  arbitrary.  Because  the  displacement  field  is  assumed 
to  vary  linearly  over  each  element  and  the  strain  displacement  relations 
are  linear,  the  strains  over  each  element  are  constant.  The  elemental 
strains  can  be  written  in  the  form: 


av^  +  CV2  +  ev^ 
bw-j  +  <i^2  9W3 

bv^  +  dv^  +  gv^  +  aw^ 
bu.|  +  du^  +  gu^ 

au^  +  CU2  +  eUj 


(3.10) 


where  Aj^  is  the  area  of  the  element,  u.,v^,w^,  (i  =  1,2,3)  are  the  u,v, 
and  w  displacements  of  the  nodes  1,  2  and  3,  respectively,  and 
a,b,c,d,e,g  are  known  constants  involving  nodal  coordinates. 

The  potential  energy  of  the  element  is  then  expressed  in  terms  of 
the  nodal  displacements  and  forces.  Minimization  with  respect  to  the 
nodal  displacements  yields  the  following  set  of  equations; 
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where  [K]  is  the  9x9  element  stiffness  given  in  Appendix  B. 

The  stiffness  matrices  of  all  the  elements  are  superposed  to  obtain 
the  global  stiffness  matrix.  Boundary  conditions  are  imposed  as 
follows  (Fig.  2): 

Displacement  Boundary  Conditions: 

V=0  along  y-0  and  W=0  along  z=0 

This  is  achieved  by  constraining  all  the  nodes  on  the  line  y=0  against 
displacement  in  the  y  direction,  and  those  on  z=0  against  displacement 
in  the  z  direction.  Due  to  the  assumed  linear  variation  of  the  displace 
ment  field,  constraining  two  adjacent  nodes  also  constrairis,  the  line  ; 
joining  them. 

Traction  Boundary  Conditions: 

T^=0  on  z=H  and  y=b 

The  traction  boundary  conditions  are  imposed  by  applying  statically 
equivalent  nodal  forces.  The  surfaces  at  z=H  and  y-b  are  stress  free 
and  equivalent  nodal  forces  are  therefore  zero. 

The  reader  is  refered  to  [1]  for  a  more  detailed  discussion  of  the 
finite  element  foimulation. 

3.3  Thermal  Loading 

The  basic  assumption  in  the  thermal  formulation  is  that  the  total 
strain  can  be  written  as  a  sum  of  a  stress  related  mechanical  strain 
and  a  free  thermal  strain. 

The  displacement  field  over  the  element  has  the  same  form  as  (3.5) 
but  the  uniform  strain  c.  instead  of  being  known,  is  treated  as  an 
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additional  unknown  that  is  common  to  all  the  elements.  It  is  equivalent 
to  the  total  laminate  strain  during  the  thermal  loading. 

The  mechanical  strain  (e)^  is  used  to  calculate  the  strain  energy 
of  the  element. 

{c}''  =  {zr  -  {c}''’  (3.12) 

where  {c}°  is  the  total  strain,  and  the  thermal  strain.  In  terms 
of  the  displacement  field  and  coefficients  of  expansion,  the  mechani¬ 
cal  strain  in  the  k^*^  layer  is: 

2  2 

^  -  (m  O.J  +  n  a2)AT 

2  2  \ 

(av^  +  0^2  ■  (f^  «V  ^  ci2)aT 

{bw.|  +  dW2  +  9W3)/A|^  -  “jAT 

(bv.|  +  dv2  +  gvj  +  +  ew2)/A|^ 

(bUi  +  du2  9^3)/^!^ 

(au^  +  CU2  +  2inn{a^  -  a2)AT 

Minimizing  the  potehtfal  energy  of  the  element  with  respect  to 
the  nodal  displacements  and  the  unknov/n  strain  c»  results  in  the 
following  set  of  equations 
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“1  ^ 
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/^x  ^ 
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^ly 

'iz 

“2 

^2x 

[K]''  ^ 

\  1 

'^2y  ^ 

^2  j 
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^2z 

"aj 

3x 

O 

1 

f ,  : 

3y 

”3 

f^ 

3z 

1  f-  -  / 

where  [K]  is  the  10x10  element  stiffness  given  in  Appendix  B. 

The  global  stiffness  matrix  is  obtained  by  the  superposition  of  the 

element  stiffness  matrices.  Boundary  conditions  are  imposed  as  in  the 

case  of  the  mechanical  load.  There  is  one  additional  equation  in  the 

thermal  problem  for  determining  the  uniform  unknown  strain  C-  It  is 

equivalent  to  the  force  equilibrium  equation  for  the  thermal  load 

n 

T.  f.  -■=  F  =  0  (3.15) 

k=l 

The  system  of  equations  is  solved  for  displacements  and,  as  in  the 
case  of  the  mechanical  loading,  the  strains  and  stresses  calculated. 

The  thermal  response  is  assumed  to  be  linear  elastic  in  this 
study.  The  stress  state  resulting  from  some  temperature  change  from 
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T.  to  is  given  by 

{a}  -f  ^  [C(T)]{^'’(T)}dT  (3.16) 

J-r  U  I 

i 

d- 

As  exact  mathematical  forms  for  [C]  and  (^1  a*”?  not  known,  continuous 

integration  cannot  be  performed,  and  for  an  incremental  solution  (3.16) 
must  be  evaluated  as  a  summation. 

nsteps 

{o}  =  Z  {Aa,}  (3.17) 

i=l  ’ 

Consider  the  ith  load  step  in  which  the  laminate  is  subject  to  a 
temperature  change  from  T.j  to  T^.  By  the  incremental  strain  theory, 
the  increment  of  stress  is  given  by 

{Aa.}|  =  [C(T)3{Ae‘'(T)}  -  (3.18) 

^  T 

However,  as  pointed  out  by  Hahn  and  Pagano  [11]  this  equation  does  not 
take  into  account  the  temperature  dependence  of  elastic  properties. 
Their  total  strain  theory  gives  the  expression  for  the  incremental 
stress  to  be 

{Aa.}l^  =  [C(T){Ac''(T)}  +  (dY  CC(T)])£‘'(T)aT  (3.19) 

The  second  term  of  (3.19)  is  difficult  to  evaluate  accurately  in  view 
of  the  limited  data  available  for  C(T).  Further,  fcr  properties  which 
do  not  exhibit  large  temperature  dependence,  the  second  term  will  be 
small  for  small  AT.  Thus  the  incremental  stresses  were  approximated  as 
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{A0.}1  :  [C(T^)]{Ae°(T^)}  (3.20) 

^2 

where  is  some  intermediate  temperature  between  T-j  and  T2,  chosen 
to  be  the  mean  in  this  study.  The  temperature  dependence  of  proper¬ 
ties  is  included  in  the  formulation  through  the  term  C(T). 

3.4  Nonlinear  Analysis 

3.4.1  Mechanical  Loading 

The  nonlinear  analysis  is  carried  out  in  an  incremental  fashion 
using  data  obtained  in  the  previous  load  step  to  calculate  the  material 
constants  for  the  current  load  step. 

Ramberg  Osgood  parameters  [23]  are  used  to  represent  the  nonlinear 
stress  strain  relations.  Typically 

n. 

e  =  |+kio’  i=l,2  (3.21) 

The  tangent  modulus  is  defined  as 


(3.24) 
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The  tangent  modulus  for  each  elastic  constant  is  calculated  assuming 
that  there  is  no  interaction  between  the  various  stresses  during  nonlinear 
behavior. 

As  indicated  in  Fig.  3  there  is  some  error  in  following  the  stress- 
strain  curve.  This  error  can  be  minimized  by  iterating  the  solution 
in  every  load  step,  or  by  using  smaller  load  steps,  as  is  done  in  this 
study. 

3.4.2  Thermal  Loading 

Temperature  dependent  elastic  properties  are  used  for  the  analysis 
of  thermal  loading.  The  elastic  moduli  E-j-j,  £22*  Poisson's 

ratios  v.|2»  V23,  etc...  and  strengths  X,  Z,  oto.,  are  input  at 
various  temperatures,  in  the  form  of  percent  retentions,  as  shown  in 
Fig.  4.  In  a  given  thermal  load  step,  the  mean  temperature  T^  is 
found  and  the  property  linearly  interpolated  between  the  nearest  higher 
and  lower  temperatures  (TpT2)  for  which  properties  have  been  input. 

These  interpolated  values  are  then  used  to  evaluate  the  stiffness 
matrices.  The  analysis  accuracy  improves  with  a  larger  number  of  input 
points,  since  the  retention  curves  for  elastic  properties  and  strengths 
are  approximated  to  greater  accuracy. 

3.5  Failure  Criterion 

The  finite  element  analysis  provides  a  three  dimensional  state  of 
stress,  presenting  a  unique  opportunity  to  study  stress  interaction 
and  failure. 

Tsai  and  Wu  [24]  proposed  that  the  failure  surface  be  represented 
in  the  form  of  a  tensor  polynomial 


20 


-  Actual  Stress  Strain  Curve 


- -  Ramberg  Osgood  approximation 

,A2  errors  in  load  steps  1 ,2 


Strain 


FIGURE  3.  DETERMINATION  OF  TANGENT  MODULI  WITH  RAMBERG- 
OSGOOD  APPROXIMATIONS 


Percent  Retention  of  Elastic  Property 
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FIGURE  4.  TYPICAL  PERCENT  RETENTION  CURVE 
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(3.25) 


The  F^j  is  a  second  order  tensor  and  F.j|^-j  a  fourth  order  tensor. 
The  numerical  values  of  the  terms  are  obtained  from  the  material 
principal  strengths.  The  tensors  simplify  greatly  for  orthotropic 
material.  The  transformations  of  the  nonzero  terms  in  these  tensors,; 
in  the  contracted  notation,  are  given  in  Appendix  C.  The  strength 
parameters  F12,  ^23  ^13  •'®Quire  special  biaxial  loading  tests 

unlike  all  other  parameters  which  can  be  obtained  from  tensile,  com¬ 
pressive  and  shear  tests.  Failure  is  predicted  to  occur  when  the  value 
of  the  polynomial  is  equal  to  or  greater  than  1.0.  The  failure  mode 
can  be  predicted  by  comparing  the  individual  contributions  of  the 
stress  components  to  the  polynomial  [25]. 


4,1  Mesh  Size 


4.  PRELIMINARY  STUDIES 

The  present  analysis  is  conducted  at  the  lamina  level,  (treating  the 
fiber  matrix  system  as  a  continuum)  and  not  at  the  micromechanical  level. 
The  finite  element  method  discretizes  the  domain  being  analyzed.  Using 
finer  grids,  one  can  get  a  better  representation  of  gradients  and  hence 
better  results.  The  problem  is  deciding  on  the  appropriate  size  of 
elements  for  the  problem  being  studied. 

Lamination  theory  results  are  accurate  in  the  interior  of  the 
laminate.  The  elements  in  that  region  can  be  much  larger  than  those 
near  the  free  edge  where  large  stress  gradients  exist  and  a  finer  mesh 
is  necessary.  The  effect  of  mesh  size  was  studied  by  using  various 
meshes  for  a  [90/0]^  laminate  that  was  loaded  with  the  same  strain  of 
0.1  percent,  keeping  all  the  material  properties  constant.  It  was 
observed  that  not  only  do  the  stresses  in  the  region  near  the  free 
edge  change,  but  the  maximum  value  of  the  tensor  polynomial  used  to 
predict  failure  changes  with  mesh  size.  A  linear  elastic  analysis 
also  predicts  different  first  failure  location,  for  the  same- laminate 
and  the  same  loading,  depending  on  the  mesh  used.  Table  1  shows  the 
location  and  maximum  value  of  the  tensor  polynomial  for  various  meshes 
(Appendix  El,  E2,  E3,  E4,  E5)  for  a  tensile  load  of  0.1  percent  strain. 
The  meshes  were  generated  using  the  mesh  generation  code  devloped  by 
Bergner  et  al  [26]. 

The  stress  distribution  is  also  a  function  of  mesh.  For  example, 
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TABLE  1 

INFLUENCE  OF  MESH  SIZE  ON  FIRST  FAILURE  IN  A 
[90/0]g  LAMINATE 


Mesh 

Elements 

Failure  Location 
on  Free  Edge 

Tensor 

Polynomial 

El 

124 

center  of  top  layer 

.238776 

E2 

230 

center  of  top  layer 

.238876 

E3 

326 

near  interface  in 
top  layer 

.2564'35 

E4 

598 

near  interface  in 
top  layer 

.266231 

E5 

878 

near  interface  in 
top  layer 

.280495 

. 

Linear  Elastic;  Applied  Strain  Loading 
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for  a  [90/0]^  laminate  exhibits  singular  behavior,  with  a  large  ten¬ 
sile  value,  at  the  free  edge  [25,27].  However,  if  the  grid  used  is  not 
fine  enough,  it  is  compressive  rather  than  large  tensile  (t’g.  5). 

In  6r/E  laminates,  there  are  approximately  20-25  filaments  through 
the  thickness  in  each  ply.  Fig.  6  shows  the  smallest  elements  in  the 
grid  drawn  on  the  photograoh  of  a  typical  Gr/E  ply.  In  the  finest  mesh 
used  in  this  study,  there  are.16  elements  through  each  ply  for  four  ply 
laminates  and  32  elements  through  the  thickness  for  two  ply  laminates. 
Therefore,  per  element,  there  is  just  over  one  filament  in  the  thickness 
direction.  For  a  laminate  aspect  ratio  of  25,  the  number  of  filaments 
calculated  to  be  in  the  smallest  element  is  approximately  3.75,  assuming 
a  fiber  volume  fraction  of  0.5.  This  mesh  (Fig.  E5)  was  modified  so 
that  it  could  also  be  used  for  a  four  layered  laminate  (Fig.  7). 

4.2  Averaging  Finite  Element  Results 

The  finite  element  formulation  used  in  this  investigation  yields 
constant  values  ror  stresses  over  each  element.  Two  adjacent  elements, 
in  general,  give  different  values  for  the  stress  at  points  on  their 
common  boundary.  In  order  to  elitriinate  the  discontinuity  of  the 
stresses,  most  finite  element  analyses  use  an  averaging  technique. 

The  following  averaging  scheme  is  used  in  this  study  (Fig.  S) . 

The  interlaminar  stresses  a-,  t  ,  t  must  be  continuous  throughout 
the  laminate.  At  a  point  A,  these  stresses  are  averaged  over  the 
elements  11,  12,  13,  14.  The  laminate  stresses  may  be 

discontinuous  across  laminate  interfaces,  but  within  each  ply,  they 
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are  continuous.  At  a  point  B  they  are  averaged  over  elements  15 
and  16. 

4.3  Linear  Elastic  Analysis 

The  tensor  polynomial  predicts  failure  to  occur  when  it  attains 
the  value  1.0.  5uppose  that  for  an  applied  strain  c.  the  stress  state 
is  determined  and  the  tensor  polynomial  calculated  as 

^j'^ij  ^  ^ijkl"ij"kl  =  “  ^  ^  .  (4.1) 

Let  failure  occur  at  a  strain  of  kc,  i.e. 

ka  +  k^3  =  1  (4.2) 

This  quadratic  equation  is  solved  for  k,  and  the  strain  at  first 
failure  determined. 

The  stresses  and  tensor  polynomial  were  evaluated  for  various 
laminates  loaded  with  an  axial  strain  of  0.1  percent.  Based  on  these 
results,  'k'  was  calculated  for  the  element  which  had  the  highest 
value  of  the  tensor  polynomial  function  for  each  laminate. 

These  values  are  presented  in  Table  2,  and  were  used  to  estimate 
the  mechanical  load  for  first  failure,  and  the  number  of  load  increments 
in  the  nonlinear  analysis. 

4.4  Stress  Free  Temperature 

During  manufacture,  laminates  are  heated  to  some  maximum  elevated 
temperature;  however,  bonding  usually  takes  place  at  some  lower 
temperature.  At  that  temperature,  the  laminate  is  in  a  stress  free 


TABLE  2 


LINEAR  ELASTIC  PREDICTIONS  OF  FIRST  FAILURE 


Laminate 

Strain  at 

First  Failure* 

[0/90].. 

0.327956 

[90/0]j. 

0.375117 

[slO]^ 

0.253761 

[*-15], 

0.22304 

[»30].  1 

0.481329 

[±45]^ 

0.476636 

H 

O 

0.456682 

[±75]^ 

0.423255 

[90/0/±45]^ 

0.318849 

[±45/0/90]^ 

0.149065 

*  Mesh  E6 
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state.  This  stress  free  temperature  is  the  reference  temperature 
for  calculating  the  residual  stresses.  depends  on  the  material 
system  of  the  laminate,  and  the  cure  cycle  used.  Tsai  [6]  suggested 
that  T^  be  experimentally  determined  from  a  [±0]  unsymmetrical 
laminate  which  warps  on  cooling.  The  temperature  at  which  the 
laminate  becomes  flat  on  reheating  is  T^.  T300/5208  is  cured  at  > 

350°F,  but  the  suggested  values  for  T^  vary  widely.  Renieri  and 
Herakovich  [1]  used  a  value  of  270°F,  while  Chamis  always  uses  the 
highest  temperature  attained  in  the  cure  cycle  as  the  value  for  T^. 

A  stress  free  temperature  of  250°F  is  suggested  iti  [18,20].  Hahn  [14] 
reheated  warped  unsymmetrical  laminates,  but  found  values  of  T^  varying 
from  250  to  300  degrees.  T^  was  choosen  to  be  250®F  for  the  present 
analysis. 

4.5  Load  Steps  for  Thermal  Load 

.  w.  A  study  was  conducted  to  evaluate  the  effect  of  cooling  the  laminate 
in  different  load  steps.  A  [90/0]^  laminate  was  chosen  because  it  ex¬ 
hibits  the  maximum  mismatch  in  expansion  coefficients  and  material 
properties.  This  laminate  was  analyzed  assuming  the  cooling  from  T^  to 
room  temperature  in  1,  2,  3,  4,  5,  6,  8,  and  10  load  steps,  and  the 
resulting  distributions  of  plotted.  Typical  variation  of  the 
stress  is  presented  in  Fig.  9.  The  largest  value  of  the  tensor  poly¬ 
nomial  was  also  determined  for  these  case  studies.  The  results 
presented  in  Table  3  show  that  the  maximum  value  decreased  with 
increasing  number  of  load  steps.  The  location  of  the  largest  tensor 
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TABLE  3 

EFFECT  OF  LOAD  STEPS  ON  THE  TENSOR  POLYNOMIAL 
[90/0]^  LAMINATE 


Maximum  Value  of 
Tensor  Polynomial 


.8767522 


.7699715 


.7252781 


.6949252 


No.  of  Load  Steps 


1 


■5 

.6705358 

6 

.6497218 

8 

.6132286 

10 

.5822445 

polynomial  value  was  independent  of  the  number  of  load  steps.  As  seen 
from  Fig.  9,  the  stress  distribution  appears  to  converge  with  increasing 
number  of  load  steps.  In  this  finite  element  analysis,  the  stiffness 
matrix  must  be  recalculated  for  each  load  step.  Using  a  grid  with  896 
elements  therefore  involves  an  enormous  amount  of  expensive  computation. 
It  was  decided  to  cool  the  laminate  in  6  load  steps  of  -30'‘F  each,  a 
compromise  between  satisfactory  convergence  and  computer  cost. 


5.  STRESS  AND  FAILURE  ANALYSIS  OF  LAMINATES 

Cross-ply,  angle-ply,  and  two  quasi-isotropic  graphite/epoxy 
laminates  were  analyzed  in  this  study.  In  order  to  obtain  the  total 
stress  state  in  the  laminate,  the  process  of  cooling  to  room  temper¬ 
ature  was  modelled  in  thermal  load  steps  with  temperature  depen¬ 
dent  properties  and  the  nonlinear  analysis  of  subsequent  mechanical 
loading  was  modelled  as  a  number  of  linear  elastic  load  steps.  Stress 
distributions  were  plotted  at  the  strain  at  which  the  first  element  was 
predicted  to  fail  using  the  tensor  polynomial  failure  criterion.  Damage 
is  predicted  to  initiate  at  this  strain.  This  study  does  not  predict 
the  ultimate  failure  strain,  but  does  predict  the  mode  of  first  failure. 
Each  Toad  step  for  the  mechanical  load  was  taken  to  be  0.05  percent 
strain,  the  choice  guided  by  the  linear  elastic  predictions  for  the 
strain  at  first  failure  in  each  laminate  and  cojnputer  cost. 

5.1  Cross-Ply  Laminates 

5.1.1  Stress  Distributions 

The  mismatch  of  the  expansion  coefficients  between  adjacent  layers 
is  maximum  in  these  laminates  and  results  in  very  high  curing  stresses. 
For  the  purpose  of  comparison,  distributions  of  non-zero  laminate 
stresses  are  presented  in  Figures  10-13  for  the  following  three  cases: 

1.  residual  thermal  stresses. 

2.  nonlinear  analysis  of  mechanical  lead  at  first  failure 
(including  residual  stresses). 
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FIGURE  13.  THERMO-MECHANICAL  c,  STRESSES  IN  A  [90/0],  LAMINATE 


3.  stresses  obtained  from  a  linear  elastic  analysis  of  pure 
mechanical  load,  scaled  to  the  first  failure  strain  as 
predicted  by  the  nonlinear  analysis. 

The  laminate  stresses  and  for  the  [0/90]^  and  [90/0]^ 
laminates  are  shown  in  Figures  10-13.  As  a  result  of  cooling,  the 
laminate  shrinks  and  is  compressive  in  the  0“  layers  and  tensile  in 
the  90°  layers,  while  is  tensile  in  the  0°  layers  and  compressive 
in  the  90°  layers.  The  stress  magnitudes  are  equal  for  both  layers 
and  stacking  sequences.  The  signs  are  opposite,  thus  satisfying 
equilibrium.  With  the  application  of  an  axial  strain  load,  stress 
reversal  occurs  for  in  the  0°  layers,  but  the  90°  layers  experience 
increased  stress  magnitude  (Figs.  10,  11).  Inclusion  of  thermal 
stresses  is  shown  to  have  a  significant  influence  on  the  overall  state 
of  stress  at  first  failure.  Comparison  of  the  linear  elastic  and 
thermo-mechanical  results  for  at  first  failure  shows  that  the  [0/90]^ 
laminate  is  preferred  with  in  the  0°  layer  being  more  than  three  . 
times  that  in  the  [0/90]^  laminate  at  first  failure. 

Edge  effects  are  seen  to  be  present  for  these  laminate  stresses. 

The  axial  stress  is  higher  in  the  boundary  layer  of  the  90° 
layers  and  lower  in  the  0°  layers.  The  transverse  stress  (Oy)  decreases 
to  zero  at  the  free  edge  for  both  laminates  and  layers  as  required  by 
the  boundary  conditions.  Careful  examination  of  the  figure  indicates 
that  the  boundary  layer  for  thermal  and  themo-mechanical  loading  is 
generally  three  to  eight  times  that  for  the  linear  elastic  analysis. 

Interlaminar  normal  stress  distributions  are  presented  in 
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Figures  5  and  14  for  the  three  loading  cases.  Moment  and  force  equili¬ 
brium  of  the  free  body  diagram  in  Figure  15  requires  that 


=  0  =>y  ^  o^dy  =  0 

Thus,  the  distribution  along  the  0/90  interface  snould  be 
equivalent  to  a  pure  couple  which  balances  the  moment  due  to  Oy. 

Since  the  sign  of  o.^  changes  when  the  stacking  sequence  is  reversed, 
the  direction  of  the  couple  should  also  be  reversed.  As  indicated 
in  Figure  14,  this  condition  is  satisfied  in  principle  by  the 
distributions  for  the  two  stacking  sequences,  for  both  thermal  and 
thermo-mechanical  loading.  The  results  in  Figure  5  for  linear  elastic 
loading  of  a  [90/0]^  laminate  also  indicate  satisfaction  of  these 
equilibrium  requirements. 

As  mentioned  in  section  4.1,  the  distribution  near  the  free 
edge  is  very  dependent  on  mesh  size.  The  general  character  of  the 
distribution  is  such  that  the  equilibrium  requirement  are  not  grossly 
violated  for  any  of  the  meshes  studied  in  this  investigation,  howeve*- 
the  magnitude  at  the  free  edge  varied  from  -80  psi  for  the  coarse 

mesh  to  +145  psi  for  the  finest  mesh  for  a  linear  elastic  analysis 


and  a^ial  strain  of  0.1  percent.  These  results  confirm  those  of  Wang 
and  Dickson  [27]  that  attains  a  tensile  value  at  the  free  edge  for 
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a  [90/0]^  as  well  as  the  [0/90]^  laminate. 

The  results  in  Figure  14  were  obtained  with  the  finest  finite 
element  iviesh  from  Figure  5  (i.e.  E5).  They  show  that  in  the 
[0/90]^  laminate  is  tensile  with  singular  behavior  for  both  thermo 
and  thermo-mechanical  loading.  It  is  also  apparent  that  the  thermal 
effects  dominate  the  boundary  layer  stress  distribution  for  the  [0/90]^ 
laminate.  Reversing  the  stacking  sequence  does  not  result  in  a  mirror 
image  of  the  stress  distributions.  Figures  5  and  14  both  indicate  a 
second  reversal  of  the  gradient  of  near  the  free  edge.  The  linear 
elastic  results  (Figure  5)  predict  a  tensile  near  the  free  edge 
whereas  the  thermal  and  thermo-mechanical  results  Show  that  thermal 
effects  and  nonlinear  behavior  have  a  beneficial  effect  on  the  magnitude 

of  a  at  the  free  edge  of  a  [90/0]  laminate. 

2  ^ 

It  is  also  apparent  from  Figure  14  that  the  significance  of 

thermal  effects  is  laminate  dependent.  Thermal  effects  dominate-  the 

a  distribution  in  [0/90]  laminates,  but  mechanical  effects  are  more 
2  ^ 

dominate  in  the  [90/0]^  laminate.  Boundary  layer  width  is,  however, 
essentially  the  same  for  both  types  of  loading  and  both  laminates 
c;ctending  over  approximately  15-20  percent  of  the  laminate  width.  The 
width  of  the  boundary  layer  in  the  [90/0]^  laminate  is  essentially  the 
same  for  all  three  types  of  loading  (Figures  5  and  14). 

Through-the-thickness  variations  of  and  for  the  residual 
stress  state  are  compared  to  the  distribution  obtained  by  Wang  and 
Crossman  [16]  using  a  linear  elastic  thermal  analysis  in  Figures  16 
and  17.  Though  the  shape  of  the  stress  distributions  is  approximately 
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the  same,  there  is  significant  difference  in  the  magnitude  of  the 
stresses.  The  maximum  value  of  in  a  [0/90]^,  for  example,  is 
predicted  to  be  2.01  ksi  by  this  analysis  compared  to  a  value  of  5,4 
ksi  from  [16],  This  difference  can  be  attributed  to  the  incremental 
analysis  using  temperature  dependent  elastic  properties.  It  should 
also  be  observed  from  Fig.  16  that,  in  all  cases,  the  maximum  positive 
value  of  occurs  within  a  layer  and  not  at  the  ply  interface. 

5,1.2  Failure  Analysis 

The  curing  stresses  in  cross-ply  laminates  are  very  high.  In  a 
[0/90]g  laminate,  the  stresses  resulting  from  cooling  the  laminate 
in  six  load  steps  were  high  enough  for  the  tensor  polynomial  to  predict 
failure.  (In  fact,  cracks  are  sometimes  observed  at  the  free  edge  of 
cross-ply  laminates  [8,20].)  For  the  purpose  of  this  analysis,  the 
[0/90]j  laminate  was  cooled  in  eight  load  ^teps  in  order  to  reduce  the 
step  size  and,  thereby,  eliminate  the  predirtion  of  failure.  All  other 
laminates  were  cooled  in  six  thermal  increments.  With  the  application 
of  mechanical  load,  first  failure  was  predicted  to  occur  at  a  strain 
of  0.05  percent  in  the  [0/90]^  laminate  and  at  0.15  percent  in  the 
[90/0]^  laminate. 

The  tensor  polynomial  is  plotted  along  the  interface  and  through 
the  thickness  for  both  laminates  in  Figures  18-20.  Failure  for  both 
laminates  was  predicted  to  initiate  in  the  90“  ply  at  the  free  edge. 
Figs.  18  and  19  show  the  variation  of  the  tensor  polynomial  along  the 
interface  in  the  90“  ply,  as  determined  from  the  curing  stresses  and 
subsequent  mechanical  loading.  The  curing  stresses  are  predicted  to 
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FIGURE  19.  TENSOR  POLYNOMIAL  ALONG  THE  INTERFACE  IN  THE 
90'’  LAYER  FOR  A  [90/0]^  LAMINATE 


FIGURE  20.  TENSOR  POLYNOMIAL  THROUGH-THE-THICKNESS  AT  THE  EDGE  OF  [0/90]^  AND 
[90/0]^  LAMINATES 
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make  a  major  contribution  to  initiation  of  failure  in  regions  close 
to  the  free  edge  In  both  laminates.  Through-the-thickness  variation 
(Fig.  20)  shows  the  effect  of  curing  stresses  to  be  significant  in 
the  90"  layers  of  both  laminates.  In  the  0»  layers  the  tensor  polynomial 
has  a  negative  value,  which  is  possible  when  using  the  tensor  polynomial 
failure  criterion.  The  maximum  value  of  the  polynomial  occurs  within 
the  90“  layer  for  both  laminates,  not  at  the  0/90  interface.  Thus, 
first  failure  is  predicted  to  occur  within  the  layer  and  not  at  the 
interface. 

The  tensor  polynomial  for  the  element  which  was  first  to  fail  was 
analyzed  in  detail;  the  individual  contributions  from  each  of  the 
contributing  stresses  are  presented  in  Table  4.  The  table  shows  that 
while  0^  make  the  largest  contribution  to  the  polynomial  at  failure 
for  both  laminates,  the  contribution  from  03  is  very  significant  in 
the  [0/90]^  laminate.  It  may  be  said  that  the  [0/90]^  laminate  fails 
in  a  mixed  03-03  mode  but  the  [90/0]^  laminate  fails  primarily  due  to 
03,  i.e.  transverse  tension.  The  [0/90]^  laminate  experiences  first 
failure  at  one-third  the  failure  strain  of  the  [90/0]^.  Since  the 
[90/0]^  is  predicted  to  fail  at  a  higher  applied  strain,  it  is  prefered 
over  the  [0/90]^  for  tensile  loading. 

5.2  Angle-Ply  Laminates 
5.2.1  Stress  Distributions 

The  angle-ply  laminates  studied  were  the  [±10]^,  [il5]j,  [±30]^, 
[±45]  [i60]j  and  [^75]^.  The  thermal  mismatch  between  adjacent  plies 


53 


TABLE  4 

FAILURE  MODE  ANALYSIS  OF  CROSS-PLY  LAMINATES 


54 


is  not  as  severe  as  that  in  cross-ply  laminates  except  for  the  [±45]^ 
laminate.  Thus,  the  residual  stresses  are  in  general  lower.  It  is 
interesting  to  note  that,  in  the  material  principal  coordinates,  the 
residual  stresses  in  the  cross-ply  and  the  [±45]^  laminate  are  the 
same,  except  at  the  free  edges.  This  is  of  course  expected  because 
of  the  tensor  property  of  the  coefficient  of  thermal  expansion. 

The  highest  absolute  value  of  each  stress  component  was  normalized 
and  plotted  versus  the  ply  angle.  Figs.  21  and  22  show  the  variation 
of  the  laminate  and  the  interlaminar  stresses  for  thermal  and  mechanical 
loading,  respectively.  The  thermal  mismatch  in  angle-ply  laminates  is 
maximum  at  45®,  and  all  thermal  stresses  (Fig.  21)  attain  their  maximum 
values  at  45”,  except  for  which  attains  its  maximum  at  30°.  This 
is  because  the  stress  state  not  only  depends,  on  the  curing  strain 
(thermal  mismatch),  but  also  on  the  elastic  modulus  and  decreases 
sharply  as  the  ply  angle  increases  from  zero,  tapering  off  at  larger 
angles  [29]. 

The  maximum  value  of  the  individual  stress  components  occur  at 
different  fiber  angles  for  mechanical  loading  (Fig.  22).  The  magnitude 
of  is  large  at  low  angles,  with  its  maximum  at  0°,  while 
attains  its  maximum  at  15°.  Three  components,  o^,,  a^,  attain 

their  maximum  value  at  30°  and  t  attains  its  maximum  at  45°.  These 
results  show  that  there  are  fundamental  differences  between  thermal 
and  mechanical  loading  of  angle-ply  laminates. 

The  distribution  of  curing  stresses  is  roughly  the  same  in  all 
angle-ply  laminates,  the  difference  being  in  the  magnitudes  of 
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NORMALIZED  STRESS 


different  ply  orientations.  The  curing  stresses  in  the  [*-^5]^  are 

typical  and  are  presented.  Fig.  23  showing  the  lamina  stresses  and 

Fig.  24,  the  interlaminar  stresses.  The  stresses  e  ,  t  ,  and  t 

y 

are  seen  to  approach  zero  as  required  by  the  stress  free  boundary 
conditions.  As  in  cross-ply  laminates,  the  curing  stresses  exhibit 
edge  effects,  with  the  presence  of  a  boundary  layer  for  y/b  greater 
than  0.9.  Fig.  24  shows  that  the  boundary  layer  for  at  the  *45 
interface  is  approximately  twice  that  of  at  the  midplane. 

Though-the-thickness  variation  of  o  and  t  near  the  free  edge 
for  a  Cj:45]g"  laminate  are  compared  to  distributions  obtained  by  Wang 
and  Crossman  [15]  in  Fig.  25.  As  in  the  cross-ply  Idminates  (Fig.  17) 
the  present  solution  predicts  much  lower  stresses.  Eoth  components 
of  stress  exhibit  sharp  gradients  in  the  vicinity  of  the  interface. 
5.2.2  Failure  Analysis 

The  interlaminar  distribution  of  the  tensor  polynomial,  at  first 
failure,  as  determined  for  thermo-mechanical  loading,  is  shown  in 
Fig.  26  for  various  fiber  angles.  This  figure  demonstrates  that  the 
edge  effects  are  dominant  at  small  angles  of  orientation,  and  that  edge 
stress  concentrations  decrease  with  increasing  angle.  At  large  angles, 
failure  is  first  predicted  at  the  free  edge,  but  elements  in  the 
interior  have  large  values  for  the  tensor  polynomial,  hence  the  entire 
laminate  is  close  to  failure.  The  tensor  pol>nonial  exhibits  a  small 
negative  value  in  the  laminate  interior  for  low  fiber  angles.  This 
is  acceptable  in  the  failure  criterion,  and  signifies  that  the  region 
is  we1>  below  failure.  These  results  indicate  tnat  the  Icmirate 
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in  an  edge  mode  for  small  fiber  angles  and  in  a  laminate  mode  for  large 
fiber  angles. 

Thermal  stresses  in  angle-ply  laminates  are  an  edge  effect.  This 
is  clearly  seen  from  Fig.  27.  where  the  tensor  polynomial  has  been 
plotted  through  the  thickness  along  the  free  edge  for  several  fiber 
angles,  for  the  curing  stresses  as  well  as  for  the  stress  state  existing 
at  first  failure.  The  presence  of  the  free  edge  and  dissimilarity  of 
material  causes  additional  stress  gradients  at  the  interface.  Failure 
is  predicted  to  initiate  at  the  interface  for  low  angles,  shift  to  the 
midplane  at  45®,  and  shift  back  to  the  interface  for  angles  greater 
than  45®.  The  maximum  value  of  the  tensor  polynomial  for  thermal 
loading  occurs  at  6=45®  where  the  property  mismatch  is  largest. 

The  stress  state  of  the  element  where  first  failure  was  predicted 
was  transformed  into  the  material  coordinate  system  and  the  individual 
terms  of  the  tensor  polynomial  evaluated  and  presented  in  Table  5. 

The  tensor  polynomial  is  completely  dominated  by  for  the  10®  and 
15°  degree  laminates,  and  the  mode  of  failure  is  therefore  predicted 
to  be  transverse  shear.  With  increasing  angle,  the  contribution  of 
decreases  while  that  from  T23  and  7^^  increase  and  the  mode  of  failure 
continues  to  be  transverse  shear  up  to  e=30®.  At  45°,  the  polynomial 
is  dominated  by  the  0^  terms,  though  there  is  significant  contribution 
from  decreases  with  further  increase  in  fiber  angle.  The 

failure  mode  for  angles  equal  to  or  greater  than  45°  is  predicted  to 
be  transverse  tension. 

Failure  was  predicted  to  initiate  at  the  free  edge  for  all 


TABLE  5 


64 


CSJ 

VO 

CO 

o 

CO 

o 

CO 

CvJ  r— 

f-* 

•cr 

to 

crt 

o 

•vr 

H 

CM 

CO 

o 

lO 

O 

o 

CM 

o 

o 

\£> 

• 

• 

• 

• 

• 

U. 

fO 

vO 

CO 

CM 

o 

CO 

VO 

CJ  #— 

o 

CO 

o 

o 

H 

o 

CO 

CO 

o 

o 

o 

m 

CTk 

CO 

to 

o 

o 

o 

LO 

• 

• 

« 

« 

• 

♦ 

Ll. 

fO 

to 

to 

o 

cy» 

LU 

CVJ  CVJ 

C3V 

o 

o 

r“ 

a: 

H 

CO 

o 

CO 

f— 

ZD 

O 

CO 

o 

o 

o 

Kl- 

• 

• 

• 

• 

* 

• 

cT 

u. 

Lu 

CVJ  cn 

o 

VO 

LO 

r— 

1- 

o 

o 

o 

CM 

o 

O 

«C 

CO 

o 

o 

o 

o 

o 

oo 

o 

o 

o 

o 

o 

o 

u. 

• 

i 

• 

• 

• 

c 

t— « 

2: 

o 

z: 

CO 

cr» 

a\ 

to 

CM 

ir> 

>- 

D 

CM 

O 

r— • 

-j 

OO 

o 

CO 

o 

U. 

o 

o 

O 

o 

a. 

• 

• 

• 

• 

oi 

* 

o 

■■lilill^B 

t/) 

CVJ  CVJ 

ur> 

VO 

CO 

VO 

D 

CM 

to 

CM 

VO 

LU 

CVJ 

o 

o 

O 

CO 

CM 

'  uo 

h- 

CVJ 

o 

o 

o 

r— • 

CM 

CM 

Li- 

• 

• 

• 

• 

* 

UJ 

1 

X 

H- 

U. 

o 

CVI 

LO 

CO 

o 

o 

tn 

D  ■ 

f- 

<K^ 

C^  1 

CO 

CM 

o 

CM 

o 

^  I 

r~~ 

5: 

U. 

O 

r-* 

VO 

VO 

a: 

• 

• 

• 

• 

• 

LiJ 

1 

1 

1 

1— 

-j  1 

cvr  r- 

iO 

CT\ 

iO 

o 

o 

<  1 

D 

O 

O 

o 

o 

o 

ZD  ! 

to 

CM 

o 

o 

o 

o 

f— 

o 

o 

o 

o 

o 

o 

U- 

• 

• 

• 

• 

• 

* 

> 

Q 

zz 

ro 

CO 

CO 

CsJ 

D 

CO 

VO 

CM 

o 

o 

CM 

CM 

r~»  ' 

o 

o 

o 

O 

O 

o 

o 

Q 

o 

■■ 

HH 

• 

1 

• 

! 

a> 

wmmmm 

4-> 

iA 

to 

XA 

fl3 

r— , 

r-n 

r~i 

C 

O 

O 

UO 

•f— 

VO 

tf 

4*1 

■f  1 

«3 

U-JI 

U~l 

L—i 

Hi 

laminates  studied.  In  the  10®  and  15°  laminates  first  failure  initiated 
at  0.003  percent  strain,  in  the  45°  at  0.0045  percent  strain  and  at 
0.004  for  the  30°,  60°  and  75°  laminates.  The  strains  at  which  first 

failure  is  predicted  is  the  same  for  some  laminates  because  the  strain 

1 

was  applied  in  load  steps  of  0.05  percent.  Linear  elastic  results  in 
[25]  indicated  that  the  [±15]^  laminate  was  the  most  critical  and  that 
the  strain  to  failure  increased  with  increasing  fiber  angle. 

5.3  Quasi-Isotropic  Laminates 

5.3.1  Stress  Distributions 

The  quasi-isotropic  laminates  analyzed  were  the  [*45/0/90]^  and 
the  [90/0/±453j.  Laminate  and  interlaminar  stresses  for  residual 
thermal  and  first  failure  under  thermo-mechanical  loading  are  presented 
in  Figs.  28-33.  Axial  (o  )  and  transverse  (o  )  stresses  in  the  90° 
layer  of  both  laminates  are  shown  in  Figs.  28  and  29,  respectively. 

The  results  show  a  strong  edge  effect  in  which  is  tensile,  and 
thus  leads  to  early  transverse  tension  failure  at  the  free  edge.  The 
width  of  the  boundary  layer  for  thermo-mechanical  loading  is  signifi¬ 
cantly  larger  than  that  for  thermal  loading.  This  is  believed  to  be 
a  manifestation  of  nonlinear  material  behavior.  The  residual  stresses 
are  shown  tc  make  a  significant  contribution  to  the  stress  state  that 
exists  when  failure  initiates. 

Various  interlaminar  stresses  are  plotted  at  different  interfaces 
of  both  laminates  in  Figs.  30  and  31.  As  indicated  in  the  figure, 
the  edge  effects  extend  further  into  the  interior  for  thermo-mechanical 
loading  than  they  do  for  purely  thermal  loading  as  was  the  case  for 
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laminate  stresses.  Though  tends  to  zero  at  the  boundary, 
and  Oj,  are  singular  at  the  free  edge.  This  was  the  case  at  all  Inter¬ 
faces  except  the  90/0  interface  where  reverses  sign  from  its  large 
negative  value,  tending  to  zero  or  some  tensile  value  at  the  free 
edge.  Such  a  behavior  was  also  predicted  for  a  [90/0]^.  laminate 
(F^gs.  5  and  14). 

Through-the-thickness  distributions  of  and  are  presented 
for  thermal  and  thermo-mechanical  loading  in  Figs.  32  and  33.  The 
stress  distributions  show  the  unloading  of  the  0*  layers  with  the 
application  of  mechanical  load.  However,  this  positive  feature  of 
thermal  stress  is  offset  by  the  fact  that  the  thermal  stress  has  the 
same  sign  as  the  stress  due  to  mechanical  in  the  90®  layers  and, 
therefore,  contributes  to  early  failure  in  that  layer.  The 
distributions  show  that  the  signs  of  the  stresses  due  to  thermal  and 
thermo-mechanical  are  the  same.  The  interlaminar  nomal  stresses  are 
predominately  compressive  for  these  two  stacking  sequences  and, 
therefore,  do  not  contribute  to  delamination.  The  only  exception 
being  near  the  midplane  of  the  [i45/G/90]j  laminate.  It  would  appear 
from  these  results  that,  of  the  two,  the  [90/0/i45],  is  the  preferred 
laminate  for  tensile  loading. 

5.3.2  Failure  Analysis 

The  distribution  of  the  tensor  polynomial  in  the  90®  layer  adjacent 
to  the  0/90  interface  of  both  laminates  is  shown  in  Fig.  34  and  the 
through-the-thickness  distributions  near  the  free  edge  are  shown  in 
Fig.  35.  Both  figures'  show  results  for  residual  thermal  and  thermo- 
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mechanical  loading.  The  90°  layer  was  chosen  for  illustration  because, 
as  indicated  in  Fig.  35,  the  polynomial  attains  its  largest  value  in 
this  layer.  The  distributions  in  Fig.  34  show  that  the  polynomial  attains 
its  maximum  at  the  free  edge.  A  boundary  layer  effect  is  very  evident 
with  the  width  at  the  boundary  layer  for  thermo-mechanical  loading 
being  several  times  that  for  thenral  loading.  This  figure  also  shows 
that  themal  effects  make  a  significant  contribution  to  the  tensor 
polynomial,  and  that  the  boundary  layer  effects  are  much  stronger  in 
the  [i45/0/90lj  laminate. 

The  through-the-thickness  distributions  in  Fig.  35  indicate  that, 
for  both  laminates,  the  maximum  value  of  the  polynomial  occurs  within 
the  90°  layer  as  opposed  to  the  0/90  interface.  In  the  [90/0/ 145]^ 
laminate,  the  maxinttjm  value  is  just  below  the  midpoint  of  the  layer 
thickness.  In  the  [*45/0/90]^,  the  maximum  value  is  at  the  center 
of  the  90°  layer,  the  midplane. 

The  terms  making  significant  contributions  to  the  tensor 
polynomial  at  failure  (in  material  coordinates)  are  presented  in 
Table  6  for  the  element  that  was  the  first  to  fail  in  each  laminate. 

These  results  show  that  there  is  significant  difference  in  the  mode 
of  first  failure  of  the  two  quasi-isotropic  laminates.  The 
[i45/0/90]j  laminate  fails  in  a  mixed  02-U3  node  whereas  the 
[90/0/i45]^  laminate  fails  in  a  predominately  02  The  influence 

iyf  the  higher  interlaminar  nonnal  stress  in  the  [r45/0/903j  is 
shown  to  lead  to  failure  at  a  lower  applied  strain. 


6.  CONCLUSIONS 


The  present  analysis  has  concentrated  on  the  evaluation  of  residual 
thermal  stresses  induced  during  curing  of  composites  and  thermo- 
mechanical  stresses  due  to  combined  thermal/mechanical  loading. 
Particular  attention  was  given  to  the  influence  that  boundary  layer 
effects  near  the  free  edge  have  on  the  initiation  of  failure  in  angle- 
ply,  cross-ply  and  quasi-isotropic  laminate  as  predicted  by  the  tensor 
polynomial  failure  criterion.  The  following  conclusions  can  be  drawn 
as  a  result  of  the  study. 

1.  Mesh  si2e  has  a  significant  effect  on  the  stress  values 
obtained  from  a  finite  element  investigation  of  the 
stress  distribution  near  the  edge  of  a  finite  width 
laminated  composite. 

2.  Thermal  effects  are  significant  in  the  boundary  layer 
of  laminated  composites. 

3.  The  boundary  layer  stress  distribution  in  jQ,ross-ply 
laminates  is  very  dependent  on  the  stacking  sequence. 

For  linear  elastic  material  behavior,  both  [0/90]^ 
and  [90/0]^  laminates  exhibit  tensile  interlaminar 
normal  str-ess  at  the  free  edge. 

4.  Failure  in  a  [0/90]^  laminate  initiates  at  approximately 
one-third  the  initial  failure  strain  of  a  [90/0], 
laminate  under  tensile  load. 
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5.  Failure  of  cross-ply  laminates  initiates  at  the  free 
edge  within  the  90”  layer,  hot  at  the  0/90  interface. 

6.  Failure  of  a  [0/90]^  laminate  is  a  mixed  mode  in 

whereas  the  [90/0]^  laminate  fails  primarily 
due  to  transverse  tension  (a^)^ 

7.  The  stress  behavior  in  angle-ply  laminates  is 

fundamentally  different  for  thermal  and  mechanical 
loading.  The  laminate  is  rwst  critical  for 

thermal  loading,  but  the  [±15]^  laminate  is  most 
critical  for  tensile  loading. 

8.  Angle-ply  laminates  with  small  fiber  angles  fail  due 
to  interlaminar  shear  whereas  laminates  with  large 
fiber  angles  fail  due  to  transverse  tension. 

9.  Failure  of  angle-ply  laminates  initiates  at  the  free 
edge.  For  small  and  large  fiber  angles,  failure 
initiates  at  the  ±o  interface.  For  intermediate 
angles,  failure  initiates  at  the  midplane. 

10.  Two  modes  of  failure  are  predicted  for  angle  ply 
laminates,  an  edge  mode  for  fiber  angles  equal  to 
or  less  than  30®  and  a  laminate  mode  for  angles 
equal  to  or  greater  than  45°. 

11.  Failure  in  the  t.t45/0/903  and  C90/0/i45]  laminates 

initiates  at  the  free  edge  near  the  center  of  the 
90®  Idyer(s).  The  [.t45/0/90l^  fails  at  a  lower 
applied  strain  in  a  mixed  mode  wnereas  the 


[90/0/i45]j  laminate  falls  primarily  in 
tension  mode. 


a  transverse 
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For  a  6  rotation  about  the  3  (z)  axis  (Fig.  1),  the  constitui.iv(‘ 
relationship  becomes 

{o}  =  [C]({e}-{u}AT) 

where 
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and  the  various  matrix  and  vector  terms  as  functions  of  the  principal 
material  properties  are  given  below  (m=cos0,  n=sin0). 


C^^+2m^n^(C^2‘^2Cgg)+n‘'C22 

n^(Cii+C22-4Cgg)+(m^+n'^)C^2 


ni 

n2 

^3 

he 

^22 

^23 

''26 


2r  s-'Jr 
m  C-|2+n  C23 


==  -nin[ni  C-i-j-n  C^^-(ni'-n  ) (^■lp'*’2Cgg)3 


'22 


12  66' 


=  n^C-,  ■,+2m^n^(Ci2+2Cgg)+m‘^G 


'll 

2  2 
n  C-j  3'rni 


'22 


-mn(n^Ci-j-m^C22)'''('0^"'^^)  (^•12^^^66^ 


'33  "  ^33 
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'44 


■'45 


'55 


mn( 023-0^3) 
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APPENDIX  B 
STIFFNESS  MATRIX 

.  Equations  (B.l)  represent  the  equilibrium  equations  for  applied 
strain  loading.  Equ's  (B.2)  represent  the  equilibrium  equations  in 
average  force  loadings.  In  these  equations,  [K]  is  the  symmetric 
elemental  stiffness  matrix,  e^iS}  and  {T1  are  force  vectors  correspond¬ 
ing  to  the  applied  strain  and  temperature  change  respectively,  (F)  is 
the  vector  of  applied  forces,  and  {xl  is  the  vector  of  unknown  nodal 
displacements. 

A.  (B.l) 

(9x9)  (9x1)  (9x1)  (9x1) 

(B.2) 

(10x10)  (10x1)  (10x1)  (10x1) 

Defining  the  following  terms 

a  =  (22-23)72 
b  =  (Y3-Y2)/2 
c  =  (23-2^)72 
d  =  (Y^-Y3)72 
e  =(2^-22)72 

g  =  (Y2-Y|)72 
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=  the  area  of  element  {l) 

'k 

F  =  average  normal  force 

where  Y-j  through  and  through  are  the  coordinates  of  the  nodal 

points  of  element  «.  in  the  Y-Z  plane,  the  element  of  the  matricies  of 
Equ.  (B.l)  can  be  defined  as  follows. 

*^11  ~  ^66®  ^22  ^  ^66^ 

>^13  ■  *  '=24  ■  <^26='=  *  ^45*"’^''* 

K,4  .  Kj5  =  (C^/  -  C45<i2,/4< 

>^15  ■  (^26“  *  545“)/''*  '<26  '  <^6®'  *  C^jdsj/A' 

K,6  ”  (E26"=  *  '45‘"!>/'''  ><27  '  'Se"'  *  '45*’"'' 

K,7  .  (Cjjba  ♦  C,5ba)M*  Kjg  +  (Cjjdc  +  C^^dcj/A' 

K,J  =  (Cjjda  t  C^jbcj/A*  K^J  =  (Cj^gc  +  C^jdej/A* 

^19  ”  ^^36®^  ^ 

"33  ■  <'555^  *  '<44  '  <'^22''^  "  '44'>^>/''' 

K3^  =  (C^gae  +  C^5gb)/A"  K^g  =  (C^^ac  +  C44bd)/A*' 

K35  =  (C^gce  +  C45bd)/A^  K^g  =  (C22ae  +  C^^bg}//^y 

•^36  =  (^26^^  C^gg^i/A^  K47  =  (C^^ba  +  C23ab)/A^ 

Sz  ~  *  C^59a)/A*'  K^g  =  (C^4bc  +  C23ad)/A 
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"38  *  *  5455':>/''‘  "49  "  (^44'’®  *  ^23"5>/''' 

"39  =  <^363®  *  Sjgej/A^ 
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+  c,/)/a' 
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+ 
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(^22^® 
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where  f s  are  nodal  forces. 

For  Equ's.  (B.2)  the  previously  defined  terms  apply  plus  the  fol 
lowing  additional  terms 

•^llO  '  ^16®  *^210  "  ^16*^  '^310  "  ^16® 

*^410  ""  ^12®  *^510  ^  ^12®  *^610  ^  ^16® 

*^710  "  ^13^  SlO  "  ^13^  SlO  "  ^13^ 

'^1010  “  Sl^ 


^10  =  ^10  - 


^  ^26V  ^  heh  ^ 

^2  ^^16^x  ■*■  ^26^y  ^  ^  ^66^xy^® 


^3  ^  ^^16^x  *  ^26^y  ^36^z  ^65^^® 


^4  =  ^‘^12"I  +  ^22"^  ^23^  ^26^^® 


^5  "  ^^12^x  *  hz^y  ^  ^23^z  ^26^xy^® 
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^6  ^  ^  ^23^.z  ^ 

h  ^  ^23^y  ^33^z  *  Ss^xy^*^ 

^8  ""  ^^IS'^x  ^23^y  '*'  ^33^z  ■*■  Ss^xy^^ 

^9  ""  ^^ISS  ^  ^23^y  ^33^z  '*'  Se^xy^^ 

^10  ""  ^12^  ^13^z  ^16^'xy^^ 

where 

e^  =  (m^a^  +  n^a2)AT 
Ey  =  (n^a|  +  m^a2)AT 

4  =  a3AT 

^xy  “  2mn(a^-a2)AT 

For  moisture  analysis  the  vector  {T}  is  identical  except  a-j. 


are  replaced  by  B-j,  B2  and  ^2 
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APPENDIX  C 

TENSOR  POLYNOMIAL  FAILURE  CRITERION 


The  tensor  polynomial  failure  criterion  in  the  contracted  tensor 
notation  (for  an  orthotropic  material  in  the  principal  material 
directions)  has  the  form 


o^+F^a^+F^aj+F^  c^+F22a2 
+2F^2«^1'’2'^^^13°1“3''’^^23°2'^3  ""  ^ 


(C-1) 


where  the  F^.  and  F^^.  terms  are  as  previously  defined  in  Chapter  3. 

In  the  xyz  (laminate)  coordinate  system,  the  tensor  polynomial  failure 
criterion  transforms  (from  the  1-2  to  x-y  by  anticlockwise  rotation  of 
+e)  into 

F^a^+F^ay+F'o^+F'Oxy+Fjia2 


'■''22'^^''33'’^''44"5z■'^55^x2 

^''66^xy^2''i6'^x"xy'-2''k‘^y"xy 

■'2^13‘^xV2''bVz  "  '' 


(C.2) 


where  the  F'  terms,  as  functions  of  the  unprimed  F's  and  e,  are  as 
follows  (m  -  x:os0,  n  =  sine) 


F]  =  m^F^+n^F2 
F^  =  n^F^+m^F2 
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=  -2mn(F^-F2) 

^11  =  ^22 
^22  ”  ^^66'*'^^12^'^  ^22 


*^33 

^^44  "  '"^*'44‘^"^^55 

•"sB  “  "^•'44‘^'^^‘'55 

^16  ”  "n^'^C2(m  Fii-n  F22)”{ff'  ■"f'  )(2Fi2+Fgg)3 

F^g  =  -nin[2(n^F^^VF22)+(m^-n2)(2F^2-'''66^^ 

*^36  '  "'™^^^13"^23^ 


^45  =  "’"(''44-^55) 

F'2  =  niV{F^^+F22-Fgg)+(mV)F^2 

'^13  '  '’’^^13'*'"^*^23 


These  are  transformations  from  the  right  handed  1-2  coordinate  system 
into  another  right  hand  coordinate  system  obtained  by  ah  anticlockwise 
rotation  of  6°  about  the  3  axis.  If  a  ply  is  oriented  at  +e°  from  the 
laminate  axis,  the  F..  are  obtained  by  using  the  above  equations  with 

*J 

the  sines  and  cosines  of  -6°. 
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USERS  GUIDE  FOR  NONCOM  3 
(Available  upon  request) 
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APPENDIX  E 

FINITE  ELEMENT  MESHES 
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N0NC0M3  FLOW  CHART 
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T300/5208  GRAPHITE-EPOXY  PROPERTIES 
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TABLE  G.2  continued 
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